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Chapter 1 

Electron and hole spin dynamics and decoherence in 
quantum dots 

D. Klauser, D. V. Bulacv, W. A. Coish and Daniel Loss 

Department of Physics and Astronomy, University of Basel, 
Klingelbergstrasse 82, CH-4-056 Basel, Switzerland 

In this article we review our work on the dynamics and decoherence 
of electron and hole spins in single and double quantum dots. The first 
part, on electron spins, focuses on decoherence induced via the hyperfine 
interaction while the second part covers decoherence and relaxation of 
heavy-hole spins due to spin-orbit interaction as well as the manipulation 
of heavy-hole spin using electric dipole spin resonance. 

I. 1. Introduction 

The Loss-DiVincenzo proposal 1 to use the spin of a single electron con- 
fined in a quantum dot as a qubit for quantum computation has triggered 
significant interest in the dynamics and control of single spins in quantum 
dots. This has led to numerous exciting experimental achievements, among 
them the realization of single electrons in single dots 2 ' 3 as well as double 
dots, 4-6 the implementation of single-spin read out, 7 ' 8 the demonstration 
of the \/SWAP operation via pulsed exchange interaction 9 and the mea- 
surement of single-spin ESR. 10 For a detailed account of the progress in 
implementing the Loss-DiVincenzo proposal, see the extensive reviews in 
Refs. 11 and 12. 

On the theoretical side, one focus was, and still is, the investigation of 
the decoherence induced by the nuclei in the host material via the hyperfine 
interaction. The first part of this review article is devoted to the discussion 
of the rich spin dynamics that results from the hyperfine interaction. We 
first give an introduction to hyperfine interaction in quantum dots (Sec. 

II. 2[) . Subsequently, we discuss dynamics under the influence of hyperfine 
interaction for the case of a single spin in a single dot (Sec. 11.30 and for a 
double dot with one electron in each dot (Sec. II. 4p . To conclude the part 
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about hyperfine interaction, we discuss the idea of narrowing the nuclear 
spin state in order to increase the spin coherence time (Sec. I1.5|) . 

The second part of the article is devoted to the dynamics and the ma- 
nipulation of heavy-hole spins in quantum dots. The motivation to study 
hole spins comes from the fact that the valence band has p-symmetry and 
thus the hyperfine interaction with lattice nuclei for holes is suppressed in 
comparison to that of the conduction band electrons. As a consequence, 
the main interest for hole spin dynamics is the relaxation and decoherence 
due to spin-orbit interaction and we discuss this in Sec. 11.61 The next task 
towards using hole spins as qubits for quantum computation is of course the 
coherent manipulation of single hole spins. A potentially powerful method 
to achieve coherent manipulation of spins is electric dipole spin resonance 
(EDSR). An analysis of EDSR for heavy holes in quantum dots will be 
presented in Sec. 11.71 

1.2. Hyperfine interaction for electrons in quantum dots 

In this part of the article concerning electron spin decoherence we assume 
that the orbital level spacing is much larger than the typical energy scale of 
the hyperfine interaction. This is the case in typical lateral quantum dots 
containing single electrons and allows one to write an effective hyperfine 
Hamiltonian for a single electron confined to the quantum-dot orbital 
ground state ipo 



where S is the spin- 1/2 operator for a single electron and 1^ is the spin 
operator for the nuclear spin at lattice site k, while v is the volume of the 
crystal unit cell and A is the hyperfine coupling strength. For GaAs, which 
is mostly used for the fabrication of lateral dots, the average hyperfine 
coupling strength weighted by the natural abundance of each isotope is 
A k, 90/j.eV. 13 In Fig. 11.11 the hyperfine coupling of the electron spin in 
a lateral double quantum dot is illustrated. The electron spin dynamics 
under TJhf have been studied under various approximations and in different 
parameter regimes. For an extensive overview, see reviews in Refs. 11, 14 
and 15. Here, we briefly mention parts of this study before we focus on a 
few cases of special interest. The first analysis of electron spin dynamics 
under f/hf in this context showed that the long-time longitudinal spin-flip 
probability is ~ l/p 2 N, 16 i.e., this probability is suppressed in the limit 
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Fig. 1.1. A double quantum dot. Top-gates are set to a voltage configuration that con- 
fines the electrons in the two-dimensional electron gas (green) to quantum dots (yellow). 
The blue line indicates the envelope wave function of the electron (blue arrow). The 
hyperfine interaction with a particular nuclear spins (red arrows) is proportional to the 
envelope wave function squared at the position of the nuclei. Thus the nuclear spins in 
the center are drawn bigger since they couple stronger to the electron spin. 

of large nuclear spin polarization p and large number N of nuclear spins 
in the dot. An exact solution for the case of full polarization (p = 1) 
gives, for both transverse and longitudinal electron spin components, a 
long-time power law decay ~ l/£ 3 / 2 by a fraction ~ 1/N on a timescale of 
r ~ hN/A ~ lfj,s (for a GaAs dot with N ~ 10 5 ). 17 The fact that this exact 
solution shows a non-exponential decay demonstrates the non-Markovian 
behavior of the nuclear spin bath. For non-fully polarized systems p < 1 
and in the limit of large magnetic fields (or high polarization p ^> l/\/~N), 
the transverse electron spin undergoes a Gaussian decay 17 ' 18 on a timescale 
r - h^N/Ay/l -p 2 (r - 5ns for GaAs with p« 1 and N ~ 10 5 ). 19 This 
Gaussian decay can be reversed using a spin-echo sequence or by preparing 
the nuclear spin system in an eigenstate of h z . 19 

Several methods to prepare the nuclear spin system have recently been 
suggested 20-22 and we discuss one of these methods 21 in Sec. 11.51 Once 
the nuclear spin system is prepared in an eigenstate, the electron spin co- 
herence is on one hand limited by dynamics in the nuclear spin system 
driven by the dipole-dipole interaction for which a worst case estimate 21 
gives ~ 100/is and on the other hand, even for an eigenstate of h z , there 
is decoherence due to the flip-flop dynamics which can be important at 
times ~ HN/A ~ 1/xs (or less, depending on the size of the electron Zee- 
man splitting). For the decay of nuclear spin polarization experiments 
suggest timescales up to tens of seconds 23-26 and hysteretic behavior of 
the nuclear spin polarization with respect to the external magnetic field 
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has been observed. 27 Further, measurements of the transport current in 
the so-called spin-blockade regime 28 revealed hysteretic behavior with re- 
spect to the magnetic field 23,29 and bistable behavior in time, 23 which is 
attributed to a bistability in the nuclear spin polarization. Further, very 
recent experiments suggest a strong dependence of the nuclear field corre- 
lation time depending on whether an electron is present in the dot or not 
and thus hyperfine mediated nuclear spin flips are a possible mechanism for 
nuclear spin diffusion. 24 This last mechanism has been estimated to lead 
to fluctuations of nuclear spin polarization on a timescale of ~ 100/zs. 21 

1.3. Single-electron spin decoherence 

In this section we look in more detail at hyperfine-induced decoherence for 
a single spin in a quantum dot in the regime of large Zeeman splitting 
e z = gusbz (due to an externally applied magnetic field b z ). If e z is much 
larger than a = [Tr{pi(h z — ho) 2 }} 1 / 2 , with ho = Tr{pih z } (where pi is the 
density matrix of the nuclear spin system), we may neglect the transverse 
term S± ■ h±_ and find that the Hamiltonian is simply 

H = (e z + h z )S z . (1.2) 

This Hamiltonian just induces precession around the z-axis with a frequency 
that is determined by the eigenvalue h z of h z , where h z \n) = h z \n). For a 
large but finite number of nuclear spins (N ~ 10 5 for lateral GaAs dots) the 
eigenvalues ft™ are Gaussian distributed (due to the central limit theorem) 
with mean ho and variance a w hA/^/N. 19 Calculating the dynamics under 
Ho (which is valid up to a timescale of ~ e z /<r 2 ~ lps, where the transverse 
terms become relevant) leads to a Gaussian decay of the transverse electron 
spin state |+) = (| T) + I l»/\/2: 19 

C° ++ (t) = -L- r dh n z e(~^ l ~)\(n\ ® (+| e (-^ot)| +) g )n) |2 

V2lT<T J-oo 

1 1 (-J?) „ . . 1 / N 2h , . 

= 2 + 2 ei ^ ) ^[(e z + ho)t]; T =- = ^ I _ p _. (1.3) 

Here again, p denotes the polarization and for an unpolarized GaAs quan- 
tum dot with N ~ 10 5 we find r ~ 5ns. Applying an additional ac driving 
field with amplitude b along the x-direction leads to single-spin ESR. As- 
suming again that e z » a, we have the Hamiltonian 



-ffESR = H + 6cos(wt)5 ; 



(1.4) 
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In a rotating-wave approximation (which is valid for (b/e z ) 2 <C 1) the decay 
of the driven Rabi oscillations is given by 30 

c m t >~ i - c+ i/5«-(5 t+ S) +o ( i 50' M 

for t > max (l/cr, 1/6, 6/2ct 2 ) and + /i - cj = 0. Here, C T E T SR (i) is defined 
in the same way as C_?_i_(t) in Eq. (|1.3[) . The time-independent constant 
is given by C = exp(6 2 /8cr 2 )erfc(6/v / 8c r ) V27rb/8cr. The two interesting 
features of the decay are the slow (~ 1/V*) power law and the universal 
phase shift of 7r/4. The fact that the power law already becomes valid after 
a short time t ~ 15ns (for b ~ a) preserves the coherence over a long 
time, which makes the Rabi oscillations visible even when the Rabi period 
is much longer than the timcscalc r ~ 15ns for transverse spin decay. Both 
the universal phase shift and the non-exponential decay have recently been 
observed in experiment. 30 In order to take corrections due to the transverse 
terms S± ■ h±_ into account, a more elaborate calculation is required. The 
Hamiltonian with flip-flop terms (but without a driving field) takes the 
form 

H B = H + ^(S+h_+S-h + ). (1.6) 

In Ref. 19 a systematic calculation taking into account these so-called flip- 
flop terms was performed using a generalized master equation, valid in the 
limit of large magnetic field or large polarization. This calculation shows 
that even for an eigenstate of h z , for which the Gaussian decay in Eq. (|1.3p 
vanishes, the electron spin undergoes nontrivial non-Mar kovian decay on a 
timescale HN/A ~ 10/is. 

Other calculations 31-33 give microsecond timescales for the electron spin 
decoherence due to electron-nuclear spin flip-flops processes. The results 
in Ref. 31 suggest that also the decoherence due to dynamics in the nu- 
clear spin system via electron mediated nuclear dipole-dipole interaction 
is suppressed by a spin echo and thus that the spin-echo decay time may 
be considerably different from the (not ensemble averaged) free-induction 
decay. 

1.4. Singlet-triplet decoherence in a double quantum qot 

We now move on to discuss hyperfine induced decoherence in a double 
quantum dot. The effective Hamiltonian in the subspace of one electron 
on each dot is best written in terms of the sum and difference of electron 
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spin and collective nuclear spin operators: S = Si + S 2 , SS = Si — S 2 and 
h= i(hi+h 2 ),5h= I (hi -h 2 ): 

H dd (t) = e z S z + h • S + 5h-SS + ^S • S - J. (1.7) 

Here, J is the Hcisenberg exchange coupling between the two electron spins. 
Similar to the single-dot case, we assume that the Zeeman splitting is much 
larger than (Sh) ms and (hi) rms , where (C) rms = [Tr{pi{0 - (0}) 2 }] 1/2 is 
the root-mean-square expectation value of the operator O with respect to 
the nuclear spin state pi. Under these conditions the relevant spin Hamilto- 
nian becomes block diagonal with blocks labeled by the total electron spin 
projection along the magnetic field S z . In the subspace of S z = (singlet 
\S), and triplet |T )) the Hamiltonian can be written as 21,34 

H sz0 (t) = ^S-S + (6h z + 6b z )6S z (1.8) 

Here, 5b z is the inhomogencity of the externally applied classical static 
magnetic field with 8b z <C e z , while the nuclear difference field 5h z is 
Gaussian distributed, as was h z in the single dot case. A full account 
of the rich pseudo-spin dynamics under H sz0 (t) can be found in Refs. 21 
and 34. Here we only discuss the most prominent features for C^g(t), 
which gives the probability to find the singlet IS*), if the system was ini- 
tialized to \S) at t = 0. The parameters that determine the dynamics 
are the exchange coupling J, the expectation value of the total differ- 
ence field xq = Sb z + 5ho and the width of the difference field as (with 
5ho = (ipi\5h z \ipi) and as = (ipi\(Sh z — fj/io) 2 ^/) 1 ^ 2 )- For the asymptotics 
one finds that the singlet probability does not decay to zero, but goes to a fi- 
nite, parameter-dependent value. 34 In the case of strong exchange coupling 
| J| » max(|xo|, as) the singlet only decays by a small fraction quadratic in 
as I 'J or x /J: 

C-(t-oo)^/ 1 - 2 ^) 2 ' \ J \»^»\^ (1 . 9) 

At short times C^g(t) undergoes a Gaussian decay on a timescale 
\/ J 2 + 4x^/4\x \as while at long times t » | J|/4ct 2 we have a power law 
decay 

CS§M ~ CS§{t - oo) + e ^ T' J '!l,r . (1.10) 
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As in the case of single-spin ESR, we again have a power-law decay, now 
with l/t 3//2 and a universal phase shift, in this case: 3n/4. Measurements 35 
of the correlator Cg z g (t) confirmed the parameter dependence of the satura- 
tion value and were consistent with the theoretical predictions concerning 
the decay. Using the same methods, one may also look at transverse cor- 
relators in the S z — subspace and find again power-law decays and a 
universal phase shift, albeit, with different decay power and different value 
of the universal phase shift. 21 Looking at the short-time behavior of the 
transverse correlators also allows one to analyze the fidelity of the ySWAP 
gate. 21 

1.5. Nuclear spin state narrowing 

The idea to prepare the nuclear spin system in order to prolong the electron 
spin coherence was put forward in Ref. 34. Specific methods for nuclear 
spin state narrowing have been described in Ref. 21 in the context of 
a double dot with oscillating exchange interaction, in Ref. 22 for phase- 
estimation of a single (undriven) spin in a single dot and in an optical setup 
in Ref. 20. Here, we discuss narrowing for the case of a driven single spin 
in a single dot, for which the details are very similar to the treatment in 
Ref. 21. The general idea behind state narrowing is that the evolution of 
the electron spin system depends on the value of the nuclear field since the 
effective Zeeman splitting is given by e z + h™. This leads to a nuclear field 
dependent resonance condition e z + /i™ — u> — for ESR and thus measuring 
the evolution of the electron spin system determines h z n and thus the nuclear 
spin state. 

We start from the Hamiltonian for single-spin ESR as given in Eq. (|1.4|) . 
The electron spin is initialized to the | j) state at time t = and evolves 
under H csr up to a measurement performed at time t m . The probability 
to find | i) for a given eigenvalue /i™ of the nuclear field operator (h z \n) = 
h"\n}) is then given by 



where S n = e z + h™ — ui and b is the amplitude of the driving field. As 
mentioned above, in equilibrium we have a Gaussian distribution for the 
eigenvalues h™, i.e., for the diagonal elements of the nuclear spin density 



matrix pi(h%,0) = (n\pi\n) = exp (-{h n z - h ) 2 /2a 2 ) j^Fhto. Thus, av- 
eraged over the nuclear distribution we have the probability Pj (t) to find 




(1.11) 



February 1, 2008 10:37 World Scientific Review Volume - 9in x 6in decrev 



8 D. Klauser et al. 

the state | J.), i.e., Pj.(i) = / dh"pi(h^,0)P"(t). After one measurement 
with outcome | J.), we thus find for the diagonal of the nuclear spin density 
matrix 36 

pi { K ,0)^p^\ h n,t m )= Pl (h n z ,0)^p^. (1.12) 

■* J. V-m) 

Assuming now that the measurement is performed in such a way that it 
gives the time averaged value (i.e., with a time resolution less than 1/6) 
we have for the probability P? of measurement result | J.) as a function of 
the nuclear field eigenvalue P™ = ^pqbrjr- Thus, by performing a mea- 
surement on the electron spin (with outcome | |)), the nuclear-spin density 
matrix is multiplied by a Lorentzian with width b centered around the /i™ 
that satisfies the resonance condition e z + ft," — uj = 0. This results in a 
narrowed nuclear spin distribution, and thus an extension of the electron 
spin coherence, if b < a . In the case of measurement outcome | \) we find 

its 1 — P n (t ) 

Pl (h:,o) ^ P ^\K,t m ) = P jW,o) 1 _ (L13) 

i.e., the Gaussian nuclear spin distribution is multiplied by one minus a 
Lorentzian, thus reducing the probability for the nuclear field to have a 
value matching the resonance condition e z + h™ — u> — 0. Due to the 
slow dynamics of the nuclear spin system (see discussion at the end of 
Sec. I1.2j) . many such measurements of the electron spin are possible (with 
re- initialization of the electron spin between measurements). Under the 
assumption of a static nuclear field during M such initialization and mea- 
surement cycles we find 

Pl (h2,o) — > P {M ^\K) = ± P i(K, o)(Pfr(i - P?) M - ai , (i.i4) 

where is the number of times the measurement outcome was | J,). The 
simplest way to narrow is to perform single measurements with b <ti a. If the 
outcome is | |), narrowing has been achieved. Otherwise, the nuclear system 
should be allowed to re-equilibrate before the next measurement. 37 In order 
to achieve a systematic narrowing, one can envision adapting the driving 
frequency (and thus the resonance condition) depending on the outcome of 
the previous measurements. Such an adaptive scheme is described in detail 
in Refs. 20 and 21. With this we conclude the part on hyperfine-induced 
decoherence of electron spins in quantum dots and move on to the heavy 
holes. 
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1.6. Spin decoherence and relaxation for heavy holes 

Now we consider the spin coherence of heavy holes in quantum dots. The 
contact hyperfine interaction between lattice nuclei and heavy-hole spin is 
much weaker than that for electrons, since the valence band has p sym- 
metry. Thus (neglecting sp hybridization) only the weaker anisotropic hy- 
perfine interaction is present. Therefore, the decoherence due to hyperfine 
interaction is suppressed for heavy holes and in this section we focus only 
on the spin decoherence due to spin-orbit interaction induced by heavy-hole 
- phonon coupling. 

From the two-band Kane model, the Hamiltonian for the valence band 
of III-V semiconductors is given by 

H hulk = H LK +i 1 J-n + H z , (1.15) 

where Hlk is the Luttinger-Kohn Hamiltonian. 38 The second term is the 
Dresselhaus spin-orbit coupling (due to bulk inversion asymmetry) for the 
valence band, 39,40 J = (J x , J y , J z ) are 4x4 matrices corresponding to spin 
3/2, fl z = P Z (P X — Py), and Q, x , Q. y are given by cyclic permutations. 
The last term in Eq. (|1.15|) Hz = — 2k^bB • J — 2<7^bB ■ J is the Zeeman 
term for the valence band 41 [k and q are the Luttinger parameters 41 and 
3 = {J xi Jy i Jz))- 

We consider a [001]-grown two-dimensional system. In the case of an 
asymmetric quantum well, due to structure inversion asymmetry along 
the growth direction, there is an additional spin-orbit term, the Bychkov- 
Rashba spin-orbit term, which, in the two-band model is given by 42 
anP x E • J, where an is the Bychkov-Rashba spin-orbit coupling constant 
and E is an effective electric field along the growth direction. Due to con- 
finement along the growth direction, the valence band splits into a heavy- 
hole subband with J z = ±3/2 and a light-hole subband with J z = ±1/2 (see 
Fig. 11.21 and Ref. 40). If the splitting A of heavy- hole and light-hole sub- 
bands is large, we describe the properties of heavy- holes and light-holes sep- 
arately, using only the 2x2 submatrices for the J z = ±3/2 and J z = ±1/2 
states, respectively. The heavy-hole submatrices have the property that 
Jx — Jy — and J z = \o z - For such a system and low temperatures, only 
the lowest heavy-hole subband is significantly occupied. In this case, we 
consider heavy holes only. In the framework of perturbation theory, 39 using 
Eq. (| 1 . 1 5[) and taking into account the Zeeman energy and the Bychkov- 
Rashba spin-orbit coupling term, the effective Hamiltonian for heavy holes 
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Fig. 1.2. Band structure of a III-V semiconductor quantum well with the [001]-growth 
direction, where E g is the band gap, A is the splitting between the light- and heavy-hole 
subbands due to quantum-well confinement, and Aso i s the splitting of the valence band 
due to spin-orbit interaction. 



of a quantum dot with lateral confinement potential U(x,y) is given by 

II = —[P; - r P:) f fl.r.//) + m;V - -'/i l'i;Ih a.. (1.16) 



-^{Pl + Py) + U(x, y) + H™ - ^g ± ^ B B ± a z 



where m is the effective heavy-hole mass, P = p + |e|A(r)/c, A(r) = 
(—yB±/2, xBj_/2, yB x — xB y \ g± is the component of the g-factor tensor 
along the growth direction, and 



HH 

SO 



= iaPla + + (3P-P + P-<7+ + 1 B-Pl<j + + H.c. 



(1.17) 



is the spin-orbit coupling of heavy holes consisting of three contributions: 
the Dresselhaus term (/3), 40 the Rashba term (a), 43 and the last term 
(7) combines two effects: orbital coupling via non-diagonal elements in 
the Luttinger-Kohn Hamiltonian (oc Pj:) and magnetic coupling via non- 
diagonal elements in the Zeeman term (oc B±). This latter term repre- 
sents a new type of spin-orbit interaction which is unique for heavy holes. 44 
Here, a = 3j a R (E z ) /2m Q A, (3 = -3-f Q rj(P^} /2m A, 7 = 37 K^B/m A, 
a± = (a x ± ia y )/2, P± = P x ± iP y , B± = B x ± iB y , mo is the free electron 
mass, 70 is the Luttinger parameter 41 (E z ) is the averaged effective electric 
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field along the growth direction of a quantum dot, and A is the splitting of 
light-hole and heavy-hole subbands. The splitting between heavy-hole and 
light-hole subbands A ~ /i~ 2 , where h is the quantum-dot height. 

The spectrum of (|1.16|) for parabolic lateral confinement [U(x, y) = 
mwj(x 2 + ?/ 2 )/2] and for vanishing spin-orbit interaction (H^§ = 0) is the 
Fock-Darwin spectrum split by the Zeeman term. 45,46 From Eq. (|1.17p . it 
can be seen that H^o leads to coupling of the two lowest states |0, ±3/2) 
to the states with the opposite spin orientations and different orbital mo- 
menta \l, =p3/2). Note that the three spin-orbit terms in Eq. (| 1 . 1 7[) differ 
by symmetry in momentum space and hence mix different states resulting 
in avoided crossings of the energy levels (see inset of Fig. II. 3[) . Due to 
this spin-orbit mixing of the heavy-hole states, the transitions between the 
states |0, ±3/2) with emission or absorption of an acoustic phonon become 
possible and this is the main source of spin relaxation and decoherence for 
heavy-holes. 40 

We consider a single-particle quantum dot, in which a heavy hole can 
occupy one of the low- lying levels. In the following, we study the relaxation 
of an n-level system, the first n — 1 levels have the same spin and the n-th 
level has the opposite spin orientation. In the framework of Bloch-Redfield 
theory, 47 the Bloch equations for heavy-hole spin motion for such a system 
in the interaction representation are given by 



is the transition rate from state j to state i, St is a constant (which has 
the value of (S z ) in thermodynamic equilibrium if R(t) = 0), 



where pure dephasing (due to fluctuations along z direction) is absent in the 
spin decoherence time Ti since the spectral function is superohmic. As can 
be seen from Eq. (|1.18[) , the spin motion has a complex dependence on the 
density matrix and, in the general case, there are n—1 spin relaxation rates. 
However, in the case of low temperatures {hqs a 3> k^T), when phonon 
absorption becomes strongly suppressed, solving the master equation, we 
find that R(t) sw 0. Therefore, there is only one spin relaxation time T\. In 
this limit, the last sum in Eq. (|1.20p is negligible and the spin decoherence 
time saturates, i.e., T2 = 2T\. 





(1.20) 
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Fig. 1.3. Heavy hole spin relaxation rate 1/Tj in a GaAs quantum dot versus an applied 
perpendicular magnetic field B± (the height of the quantum dot is chosen to be h = 5nm, 
the lateral size Iq = y/ h/muiQ = 40 nm, k = 1.2, 70 = 2.5, g± = 1. 5, 49 and the 
other parameters are given in Ref. 40 ). Inset: Energy differences of lowest excited levels 
with respect to the ground state ^0^+3/2- The second avoided crossing comes from the 
spin-orbit interaction and the in-plane magnetic field B|| (3rd term in Eg. I1.17)l . The 
anticrossing gap is proportional to , implying that the coupling between corresponding 
states can be controlled externally. 



Note that in contrast to electrons 48 there are no interference effects 
between different spin-orbit coupling terms, thus the total spin relaxation 
rate 1/Ti is the sum of rates l/T x = 1/Tf + l/Tf K + i/t|: 40 ^ 44 

1 , ■ 



T 1 BR z \3lj + + ojz - loz 
— jj cx (3 u>z I 

±1 \U! + +UJz UJ-—UJZ 

— ex 7 2 g la* ± + = , 1.21 

j>ll 11 \2u) + +ujz 2u— — Uz/ 



where w± = ^/lo^ + lv 2 /A ± w c /2, uj z = g±fJ,BB±/h, B\\ = ^JB 2 + B 2 . In 
Fig. 11.31 the total spin relaxation rate 1/Ti is plotted as a function of per- 
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pcndicular magnetic field B± . There are three peaks in the relaxation rate 
curve at toz = w_ , 2lu- , and 3lu- , which are caused by strong spin mixing 
at the anticrossing points. In the inset, the first (third) avoided crossing 
resulting from Dresselhaus (Rashba) spin-orbit coupling corresponds to the 
first (third) peak of the spin relaxation curve in Fig. 11.31 At non-zero in- 
plane magnetic fields (-By), there is an additional peak which is due to an 
anticrossing between the energy levels E Q+3 / 2 anc > £^2,-3/2 ( see the second 
avoided crossing in the inset). Note that the spin relaxation rate for heavy 
holes is comparable to that for electrons 40,50 due to the fact that spin-orbit 
coupling of heavy holes is strongly suppressed for flat quantum dots (see 
Eq. (JTTT7J)), as confirmed by a recent experiment. 51 

1.7. Electric dipole spin resonance for heavy holes 

Let us now consider methods for the manipulation and detection of the 
heavy- hole spin in quantum dots. For electrons in two-dimensional struc- 
tures, an applied oscillating in-plane magnetic field couples spin-up and 
spin-down states via magnetic-dipole transitions and is commonly used in 
electron spin resonance, Rabi oscillation, and spin echo experiments. 10 It 
can be shown that magnetic-dipole transitions (An = 0, A£ — 0, and 
As = ±1) are forbidden and, due to spin-orbit mixing of the states |0, ±3/2) 
with iflf 1 1,^3/2), electric-dipole transitions (An = ±1, A£ = ±1, and 
As = 0) are most likely to occur. Therefore, the heavy holes are affected 
by the oscillating electric field component and not by the magnetic one. 

We consider a circularly polarized electric field rotating in the XY-plane 
with frequency to: E(£) = _E(sinu;i, — coswt, 0). Therefore, the interac- 
tion of heavy holes with the electric field is described by the Hamiltonian 
H E (t) = (\e\E /mui)(cosujtP x + s'mu}tP y ). The coupling between the states 
|±) is given by (+\H E {t)\-) = H E _ = (H E + )* = d so Ee- luJt , where 

d so = (\e\l/2u)(fcu- + 0+u + ) (1.22) 

is an effective dipole moment of a heavy hole depending on Dresselhaus 
spin-orbit coupling constants, perpendicular magnetic field B±, lateral size 
of a quantum dot, and frequency uj of an rf electric field. 

In the framework of the Bloch-Rcdficld theory 47 (taking into account 
also off-diagonal matrix elements) , the effective master equation for the den- 
sity matrix p nm assumes the form of Bloch equations, 47 with the detuning 
of the rf field given by <5 r f = lu% — uj. 2dsoE/H is the Larmor frequency, 
T\ = 1/(W_| + W [_) the spin relaxation time (W nm is the transition 
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rate from state m to state n), T 2 — 2Ti 40 the spin decoherence time, and 

p^ = (W-\ — W-+)Ti the equilibrium value of p z without rf field. 

The coupling energy between a heavy hole and an oscillating electric 
field is given by 



(H E (t)) =Tr( /9 7J £ (t)) = -d so -E(i), 



(1.23) 



where dso = dso(ip h — ip~\ , P-\ + P h 0) 1S the dipole moment of a 

heavy hole. Therefore, the rf power P = —d(H E (t))/dt = —ujdsoEp- 
absorbed by a heavy-hole spin system in the stationary state is given by 52 

2oj(d so EfT 2 p T Jh 
l + 8 2 s T 2 + (2d so E/h) 2 T 1 T 2 1 • ' 

In Fig. 11.41 the dependence of P on a perpendicular magnetic field B± 
and frequency w of the oscillating electric field is plotted. The rf power P 
absorbed by the system has three resonances and one resonant dip. The 
first resonance appears when the energy of rf radiation equals the Zeeman 
energy of heavy holes: B^ 1 — huj/g±pB- The shape of this resonance (at 
certain uj) is given by P w hup T z /2h[l + h 2 S^T 2 / {2d so E) 2 T 1 }. 
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Fig. 1.4. Absorbed power P (meV/s) as a function of perpendicular magnetic field B± 
and rf frequency u) (T2 = 2Ti, E = 2.5 V/cm, Bu = 1 T, and the other parameters are 
the same as those in Fig. 11.31 1 . 



If the first and second resonances are well separated (uj -C then 
the absorbed power can be estimated as 

P » 2uj(d S oE) 2 p]:/hS 2 { T 2 (1.25) 
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in the region of the second and third resonances and the resonant dip. The 
second resonance corresponds to an anticrossing of the levels -Eo,-3/2 aid 
£i, 3/2 (see the first avoided crossing in the inset of Fig. II. 3| at w_ = luz 40 
(at B = B 1 ^ ). At the anticrossing point, there is strong mixing of the 
spin-up and spin-down states and the dipole moment of a heavy-hole 
spin system is maximal rf™Q X = \e\luiz/2u} and equals half of the low- 
est electric dipole moment of a quantum dot (\e\loj ^ / lo) . Therefore, the 
height of the second resonance is given by (elu>zE) 2 /2hcu6^ f T2. The res- 
onant dip appears at B^_ = (Hujq /2g±/j,B ) \J 2mo /g±m, which corresponds 
to /3-f w_ + Pt UJ + = an d to zero dipole moment (see Eq. (|1.22j) ). The 
third resonance reflects the peak in the spin decoherence rate T 2 _1 due to 
an applied in-plane magnetic field (sec Fig. II. 3[) at the second anticross- 
ing point (the second avoided crossing in inset of Fig. 1 1 . 3[) at 2lj_ = ojz 
(B 1 ^ 3 = Ahioo/ 9±^B\/^ + 4mo/fli"j)- From the positions of the resonances 
we can determine g±, m, and u>o, from the shape and the height of those 
we can extract information about the spin-orbit interaction constants a, (3, 
and spin-orbit interaction strength due to in-plane magnetic field (which 
is proportional to 70 k/ A). Moreover, we can determine the dependence of 
the spin relaxation and decoherence times on B± . 



1.8. Conclusions 

We have discussed the rich dynamics of single electron spins in single and 
double quantum dots due to hyperfine interaction with the nuclei. Key fea- 
tures are non-exponential decays of various kinds and a remarkable univer- 
sal phase-shift. Further, we have studied spin decoherence and relaxation of 
heavy holes in quantum dots due to spin-orbit coupling. The spin relaxation 
time T\ for heavy holes in flat quantum dots can be comparable to that for 
electrons 40,50 as confirmed by experiment. 51 The spin decoherence time for 
heavy holes is given by T2 = 2T\ at low temperatures. There is strong spin 
mixing at energy-level crossings resulting in a non-monotonic dependence 
T\(B). We have proposed a new method for manipulation of a heavy-hole 
spin in a quantum dot via rf electric fields. This method can be used for 
detection of heavy- hole spin resonance signals, for spin manipulation, and 
for determining important parameters of heavy holes. 
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